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^ ABSTRACT 

O 

^-H We show that the 4.44 and 6.13 MeV line emission observed with COMPTEL from 

■^ Orion is consistent with gamma ray spectra consisting of a mixture of narrow and broad 

^ lines or spectra containing only broad lines. We employed several accelerated particle 

'q^ compositions and showed that the current COMPTEL data in the 3-7 MeV region alone 

Q cannot distinguish between the various possibilities. However, the COMPTEL upper limits 

.<— > in the 1-3 MeV band favor a composition similar to that of the winds of Wolf-Rayet stars 

Cd of spectral type WC. The power dissipated by the accelerated particles at Orion is about 
4 X 10 erg s~ . These particles are not expected to produce significant amounts of Al. 



Submitted: August 10, 1994 
Accepted: October 17, 1994 



Subject headings: gamma rays: theory — Galaxy: abundances — stars: Wolf-Rayet- 
nuclear reactions, nucleosynthesis, abundances — 
acceleration of particles — stars: sup ernovae: general 



In press. The Astrophys. J. Letters 



1. INTRODUCTION 

Gamma ray line emission in the 3 to 7 MeV range was observed from the Orion complex 
with COMPTEL on the Compton Gamma Ray Observatory (Bloemen et al. 1994). The 
radiation shows emission peaks near 4.4 and 6.1 MeV, consistent with the deexcitation of 
excited states in C and O. These lines can only be produced by accelerated particle 
interactions (Ramaty, Kozlovsky and Lingenfelter 1979). Even though gamma ray lines 
produced by accelerated particles have been seen from many solar flares (e.g. Chupp 1990), 
this is the first instance that such lines are observed from a galactic source. The Orion 
spectrum, however, is very different from solar flare spectra. Whereas Orion reveals only 
upper limits below 3 MeV, solar flare spectra show strong emission in this region, both in 
lines (mainly due to Ne, Mg, Si, S and Fe) and electron bremsstrahlung (Murphy et al. 
1991). 

Bloemen et al. (1994) suggested that the excitations are produced by the accelerated 
C and O nuclei themselves, rather than by accelerated protons and a particles bombarding 
ambient C and O. Bykov and Bloemen (1994) proposed that the interactions are predom- 
inantly due to nuclei around 30 MeV/nucleon accelerated by colliding stellar winds and 
supernova explosions, and suggest that abundances resulting from supernova explosions 
could explain the COMPTEL upper limits below 3 MeV. Nath & Biermann (1994) also 
pointed out that the stellar winds of the young O and B stars in Orion will produce shocks 
which are efficient particle accelerators. Clayton (1994) suggested that the same particles 
which produce the gamma rays in Orion will lead to an enhanced Al production which 
could be sufficient to account for both the total galactic Al mass implied by gamma ray 
observations (e.g. Diehl et al. 1993) and the Mg excess seen in meteorites in the early 
solar system (Lee, Papanastassiou & Wasserburg 1977). 

In this letter we present results obtained by employing a detailed nuclear deexcita- 
tion code to analyze the COMPTEL Orion observations. This code, originally developed 
for application to gamma ray line production in the interstellar medium (Ramaty et al. 
1979), has been extended and used to analyze solar ffare gamma ray emission (Murphy et 
al. 1991). The code incorporates a large number of nuclear reactions, and it allows us to 
calculate deexcitation line spectra in either a thin or a thick target model, for various ambi- 
ent and accelerated particle compositions, and various accelerated particle energy spectra. 
We address the following questions: do the 3-7 MeV data really favor an accelerated par- 
ticle population at Orion that is enriched in C and O relative to protons and a particles; 
what are the implications of the absence of observed emission below 3 MeV on the particle 
composition and energy spectrum; what information do the shapes of the lines provide on 
the spectrum of the accelerated particles; what other gamma ray lines could be observed 
from Orion; what is the total power in accelerated particles; and is the Al production in 



Orion significant? We also comment on the differences between tfie solar flare and Orion 
gamma ray spectra which might yield clues on the nature of the acceleration mechanism. 

2. ANALYSIS 

To calculate the gamma ray line production we must make assumptions on both 
the properties of the accelerated particles and the interaction model. The accelerated 
particle spectra and compositions predicted by various acceleration theories depend on 
many parameters, especially in the low energy range ( ^ 100 MeV/nucleon). Since the 
observed gamma ray lines from Orion are predominantly due to particles in this energy 
range, for simplicity we have chosen a spectral form which characterizes the source spectra 
Ni{E) of the various accelerated particle species in terms of a single parameter. We take 
Ni{E) = const for E < E^ and oc E at higher energies. While other spectral forms are 
also possible, we show that by varying Ec over a broad range (2-100 MeV/nucleon), we 
are able to fully explore the dependence of main features of the gamma ray emission (i.e. 
its spectrum and normalization) on the hardness of the accelerated particle spectrum. 

The other important property of the accelerated particles is their composition which 
depends on not only the acceleration mechanism but also the ambient medium from which 
the particles are accelerated. Rather than relying on the predictions of acceleration mech- 
anisms, we take the ratios of the iVj's at the same energy per nucleon proportional to any 
one of the four compositions shown in columns 3 through 6 of Table 1. We assume that 
the composition of the ambient medium is that of the solar photosphere given in column 
2 of this table. The cosmic ray source composition (CRS) is from Mewaldt (1983) except 
H for which we took the solar value relative to He. The SNISM© and SN25M0 composi- 
tions represent the ejecta of the ISMq and 25M0 supernovae given by Weaver & Woosley 
(1993). The WC composition represents the extreme late phase winds of Wolf-Rayet (WR) 
stars of spectral type WC. In such stars He, C and O are found in great abundance but 
there is no evidence for H (Abbott & Conti 1987). The WC values given in Table 1 for 
^He, ^^C, l6o, 20^e, ^^Ne and ^SMg are based on calculations (Maeder & Meynet 1987, 
figure 13) for a 6OM0 star. The other WC abundances given in Table 1 are solar values 
normalized to Ne, except H for which we arbitrarily set the abundance equal to that 
of O. In this WC composition, accelerated H plays essentially no role in gamma ray 
production. 

We evaluate the gamma ray line production in a thick target model in which the 
accelerated particles produce gamma rays as they are slowed down by Coulomb energy 
losses in a neutral medium until their energies fall below the thresholds of the various 
nuclear reactions. This model is energetically more efficient than a thin target model in 
which the particles escape from the target region after having lost only a small fraction of 
their energy. Furthermore, because the rate of energy loss increases with increasing nuclear 



charge, the contribution of the accelerated Ne-Fe nuclei, which produce gamma ray lines 
mostly below 3 MeV, is suppressed in a thick target relative to the contributions of C and 
O which produce lines mainly above 3 MeV. 

Bloemen et al. (1993) have presented an upper limit on the ratio of the 1-3 MeV 
to 3-7 MeV emission, R = 2 J^^^^ E^Q{Ej)dEj/ J^^^^ E^Q{Ej)dEj < 0.13(2cr), where 
Q{Ej) is the differential gamma ray emissivity in Orion. We have calculated R for the four 
accelerated particle compositions and the results are shown in Fig. 1 (panel a) together 
with the COMPTEL upper limit. We see that the WC composition yields very low i?'s 
which are consistent with the data for practically any accelerated particle spectrum. This 
is because the 1-3 MeV band is populated mostly by nuclear lines from Ne-Fe, and these 
nuclei are essentially absent in the WC composition. The other compositions yield higher 
i?'s. The discrepancy between the data and both the CRS and SNISMq values is greater 
than 3(7 at all -Ec's, and increases rapidly below 5 MeV/nucl. The SN25M0 values are also 
inconsistent with the data, except for Ec^'s near 5 MeV/nucl. 

In panel (b) of Fig. 1 we show the ratio of the gamma ray line yield to the total energy 
in accelerated particles for the four compositions. Here W = T^iAi J ENi{E)dE, where Ai 
is atomic number. We see that the efficiency increases with increasing Ec and that it is 
the highest for the WC composition, which is more efficient, by a factor of about 7, than 
the other compositions. 

In Fig. 2 we compare the COMPTEL data with calculated spectra obtained by con- 
volving calculated high resolution spectra with a Gaussian (a = 0.125 MeV) that approx- 
imates the energy resolution of the COMPTEL instrument. Both the CRS and SNISM© 
spectra (panel a) provide acceptable fits to the 3-7 MeV data, although, as shown in Fig. 1, 
these spectra overproduce in the 1-3 MeV band. In panel (b), the SN25M0 spectrum is 
probably unacceptable because it does not produce enough emission around 4.4 MeV. On 
the other hand, the WC spectrum in panel (b) fits the data quite well, and, as shown in 
Fig. 1, this spectrum is also consistent with the 1-3 MeV upper limit. 

In Fig. 3 we show our calculated high resolution spectra for the CRS and WC compo- 
sitions normalized to the observed fiux in the 3-7 MeV range (Bloemen et al. 1993). The 
CRS spectrum shows a variety of narrow lines superposed on broad lines. The fact that 
this spectrum fits the 3-7 MeV data quite well (Fig. 2 panel a) shows that it is not neces- 
sary to suppress the narrow line component as suggested by Bloemen et al. (1993). The 
strongest predicted narrow line fiuxes for the CRS composition are given in Table 2, where 
the indicated energy ranges contain the predicted photon fiuxes. The WC spectrum (Fig. 3 
panel b) consists mostly of broad lines, particularly those of C and O. The strongest 
narrow line is at 1.275 MeV from Ne whose predicted fiux is about 10 photons cm 
s~ . The confirmation of the absence of narrow Ne-Fe lines would provide support for 



our suggestion that the observed broad C and O hues from Orion indeed originate in WR 
stars. 

We have also calculated the 0.511 MeV line flux resulting from positron production 
by the accelerated particle interactions (Kozlovsky et al. 1987). Taking the positronium 
fraction equal to 0.9 (Guessoum, Ramaty & Lingenfelter 1991) we obtain fluxes of 2 x 
10 and 1.5 x 10 photons cm s~ for the CRS and WC compositions, respectively. 
The 0.511 MeV line and its associated positronium continuum could probe the physical 
conditions of the annihilation site (Guessoum et al. 1991). 

The widths of the broad lines depend on Ec- We found that for the WC composition 
Ec should not exceed about 20 MeV/nucleon because otherwise the broad 4.4 MeV line 
will no longer flt the data. Such an upper limit cannot be placed on the spectra which 
also include narrow lines (e.g. the CRS spectrum). However, if Ec, exceeds about 50 
MeV/nucleon then spallation reactions will produce a feature around 5.2 MeV (from N 
and O) for which there is no evidence in the data. 

In Fig. 4 we show the ratio of the Al yield to the 3-7 MeV photon yield based 
on cross sections given by Kozlovsky et al. (1987). Clayton (1994) estimated an Al 
production in Orion of about 5 x 10 nuclei per 10 years. Using the results of Fig. 4 
for Ec >5MeV/nucl and the observed 3-7 MeV flux of 10~ photons cm~ s~ , we obtain 
production rates which are lower by factors of 10 to 300, depending on the composition. 
It is thus unlikely that the accelerated particles in Orion play a major in Al production. 

3. DISCUSSION AND SUMMARY 

We have shown that the observed Orion data in the 3-7 MeV range can be flt by 
gamma ray spectra which consist of narrow and broad lines or by spectra which contain 
only broad lines. A typical example of the former is the spectrum produced by accelerated 
particles with the cosmic ray source (CRS) composition. On the other hand, to produce a 
gamma ray spectrum consisting of only broad lines, it is necessary to suppress the protons 
and a particles relative to the accelerated heavier nuclei. Such a suppression has never 
been encountered in any of the well studied accelerated particle populations, e.g. solar 
energetic particles (SEP) or cosmic rays. Nor is such a suppression predicted by known 
acceleration theories. As we have suggested, a possible scenario is the acceleration of the 
C and O rich winds of WR stars by shocks. However, none of the known WR stars (e.g. 
Lang 1991) is in Orion. But it is possible that some WR stars are hidden in or behind 
clouds or were active there in the near past. 

We found that the compositions of the ejecta of 15M0 and 25M0 supernovae do not 
offer any advantages. The former is in conflict with the 1-3 MeV upper limit at about the 
same level as the CRS composition, while the latter, because of its low C abundance, does 
not produce a good flt to the 3-7 MeV data (see Fig. 2, panel b). 



For the WC composition with Ec = 20 MeV/nucleon, the 3-7 MeV gamma ray yield is 
about 6 photons per erg (Fig. 1 panel b). Within the present constraints this yield is close 
to the maximum. The observed 3-7 MeV flux of 10~ photons cm~ s~ then implies that 
the power dissipated by the accelerated particles at Orion is about 4 x 10 erg s~ . In 
the WC scenario, for a mass loss of 5 x 10 Mq per year and a wind velocity of 4000 km 
s~ (near the maximal values given by Abbott & Conti 1987), the available kinetic energy 
is about 2.5 x 10 erg s~ per star. If the conversion efficiency of wind kinetic energy 
into shocks and accelerated particles is close to unity, a few WR stars could account for 
the observed gamma ray line emission in steady state. But such stars are presently not 
observed. Alternatively, if the wind particles are accelerated over the entire WC phase of 
^10 yrs and a significant fraction of them subsequently diffuse into contact with dense 
clouds where their stopping time is much less than their acceleration time, then a single 
recent WR star could have produced all of the accelerated particles responsible for the 
emission and no longer be visible. 

We pointed out that the solar flare gamma ray spectra show much higher ratios 
of 1-3 MeV to 3-7 MeV fluxes than does Orion. It was shown (Murphy et al. 1991) 
that this enhanced emission below 3 MeV is, in part, due to the enrichment of the flare 
accelerated particle population in heavy nuclei. Such enrichments are routinely seen in 
direct observations of SEP from impulsive flares (e.g. Reames, Meyer & von Rosenvinge 
1994). These impulsive SEP events are also rich in relativistic electrons. On the other 
hand, in gradual SEP events the composition is coronal and the electron-to-proton ratio is 
low. The acceleration in impulsive events is thought to be due to gyroresonant interactions 
with plasma turbulence while in gradual events it is the result of shock acceleration. The 
fact that the ratio of bremsstrahlung-to-nuclear line emission in Orion is very low lends 
support to the shock acceleration scenario. 

We wish to acknowledge the role of R. J. Murphy in the development of the nuclear 
deexcitation line code and thank A.G.W. Cameron for reminding us of the contribution of 
Na to Al production. We also wish to acknowledge flnancial support from the US-Israel 
Binational Foundation and NASA under Grant NAGW 1970. 
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TABLE 2. Predicted Narrow Line Fluxes 
for the CRS Composition Normalized to the 
Observed 3-7 MeV Flux 



Energy Range 


Isotope 


Photons 


(MeV) 




cm-2 s-i 


0.40-0.50 


^Li-^Be 


1.8 X 10-^ 


0.475-0.48 


^Be 


3.1 X 10-6 


0.84-0.854 


56Fe 


5.5 X 10-6 


1.232-1.244 


56Fe 


3.4 X 10-6 


1.61-1.66 


20Ne 


6.7 X 10-6 


4.35-4.52 


12 Q 


1.6 X 10-^ 


6.0-6.25 


160 


1.2 X 10-^ 



0.6 
0.5 

^ 0.4 

I 

d 0.3 

I 

^0.2 
0.1 

0.0 





1 L 1 1 • 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ' ' 1 1 1 1 1 1 1 1 1 

\sN15Mo 

\ Thick Targ 
\ Ambient Medium 

\ \ (a) 


1 ' 1 ' 1 1 1 II 

et 

: Solar 

CRS 


XcRS 






^^\ ^^^ 


SN25Mo_ 

1,1, 


^^ COMPTEL 
V 

1 ,1.1 1 1 1 1 1 , 1 . 1 .1,1,1. 



00 



101 



4 5 



102 



0) 



2 

01 



m 2 

o 100 

+-> 

o 

^ 2 



^ 10-2 



CO 10-3 



-] — 1 — I — ■— r 



T 1 — I — I — r 



I ' I ' I ' I ' I 1 — I — err 



Thicl^ Tar 



100 




CRS 
C 

■ ■ ■ ■ SNlSMo 
SN25Mo 

Ambient Medium: Solar 



3 4 5 101 2 3 4 5 

Ec(MeV/nucl) 



102 



Fig. 1. (a): Ratio of 1-3 to 3-7 MeV emission for various compositions; R and Ec are 
defined in tfie text, (b) Tfie 3-7 MeV gamma ray yield per erg deposited for the various 
compositions. 
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Fig. 2. Fits to the COMPTEL data for the various compositions. 
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Fig. 3. Calculated high resolution spectra. The comparison of the WC (Wolf-Rayet) and 
CRS (cosmic ray source) cases shows the effects of the suppression of protons and Ne-Fe 
relative to C and O in the accelerated particle population. 
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